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ABSTRACT 


Two new approximate algorithms with O(klogk) for 
the rectilenear Steiner tree are proposed. Both algo- 
rithms base upon the method which makes minimum 
spanning tree on the modified Delaunay net with the 
triangular Steiner points as the more virtual generat- 
ing points, because each point and each edge weight 


_ on the modified Delaunay net have been fixed. 


1 INTRODUCTION 


The rectilinear Steiner tree problem (RSTP) is the 
‘Steiner tree problem in the plane where instead of the 
familiar Euclidian metric L2,the distance between any 

pair of points (21,41) and (x2, y2) is given by the rec- 


tilinear metric Ly (i.e.,] 21 — v2 | +| 91 — ye |). 
The Steiner tree problem in the Euclidean plane 


_ (ESTP) is defined : find the shortest network span- 
ning a set of given points in the Euclidean plane. In 
order to obtain such a network, it is usually necessary 


to introduce additional points, called Steiner points, 


at which the edges meet one another. 
_ Among a variety of practical applications of RSTP, 
- one could mention wire layout for circuit boards, and 


layout of electrical and mechanical system that in 
building or along city streets[8]. 

However, it has shown that this problem (RSTP) is 
also NP-complete[1]. That is, it has became necessary 
to obtain polynomial-time approximation algorithms 
and, as a consequence, may approximation algorithms 
have been proposed{2]-[8]. . 

In Ref.[6],the new algorithm with O(k*) is proposed 
by using the structural propertree of the paths in the 
RST and the concept of triangular Steiner points on 
Delaunay net. | 

Meanwhile, Refs.[4Jand[5] gives the algorithms with 
O(klogk) for RSTP. In this paper, we propose two new 
algorithms with also O(klogk) one of which is, so to 


speak, the extended version of Ref.[6]. 


2 PRELIMINARIES 


Let G = (V, £) be an grid graph. Then V denotes the 
set of nodes of G, with | V |= n, and E denotes the 
set of edges of G. In this paper we assume that the 
weight of each edge of the grid graph G is equal to 1. 
Let P be a finite set of generating points in the grid 
graph such that | P |= k. Let T(V', E’) or T be an 
RST,where V’ C V and E'C E. 

For the Steiner point in the grid graph, the following 
two properties have been shown. 

[Lemma 1][3] 

All Steiner points either have the degree three or the 
degree four. a 

[Lemma 2][2] 

The maximum number of Steiner points in T is at 
most k — 2. E 

One of the most useful data structures in computa- 
tional geometry is the Voronoi diagram. 

(Definition 1)[2] | 

The Delaunay net is defined as a dual graph of the 
Voronoi diagram. Then the Delaunay net is repre- 
sented by the set of triangles whose apexes are the 
given generating points, and this fact is called the De- 
launay triangulation. | 


3 TRIANGULAR 
STEINER POINT 


We consider the optimum Steiner tree with three gen- 
erating points (this tree is denoted by the 3-point op- 
timum Steiner tree or T3 hereafter) in the grid graph 
G and we consider a triangle of Delaunay net. Since 
the number of the apexes of a triangle is three and 
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the three sides of a triangle are corresponding to the 
edges of Delaunay net, we can apply the properties of 
the 3-point optimum Steiner tree to each triangle of 
the Delaunay net. First, for the Steiner point in the 
3-point optimum Steiner tree, the following lemma can 
be obtained from Lemma 2. 
[Lemma 3] 
The optimum Steiner tree with three generating points 
has at most one Steiner point on the grid graph. | 
Here, we define the triangular Steiner point. 
(Definition 2) 
The triangular Steiner point is defined as the Steiner 
point of the optimum Steiner tree whose generating 
points are the three apexes of each triangle in the De- 
launay net. That is, consider the triangle with the 
three points P;(xi, yi), P;(wj,y;),and P(r, yx), where 
vj; < ui < vy and yi < ye < Yj. 
Then the triangular Steiner point is given as 
Ps(xi, yk). ‘ 
Let S be the set of the triangular Steiner points. i 
_ Furthermore, we define the modified Delaunay net. 
(Definition 1’) 
The modified Delaunay net is defined as the Delaunay 
net whose vertexes set is PUS. i 


4 APPROXIMATION AL- 
GORITHMS WITH O(k 
log k) 


[MAIN ALGORITHM] 
Step A) Make a Delaunay net from P; 
/* P ; the set of generating points */ [O(klogk)] 
Step B) Make triangular Steiner points of the Delau- 
nay launay triangulation; . 
/* S ; the set of triangular Steiner points */ [O(k)] 
Step C) Make a modified Delaunay net from U ; 
/*U;U =PUS */ [O(klogk)] . 
Step D) Make an MST of the modified Delaunay 
net by using [ALGORITHM 1] or [ALGORITHM 2] 
[O(klogk)] 
Step E) Delete the triangular Steiner points of 
degree(Si) < 2; /* Si; Si € 5 */ [O(K)] 
Step F) Move all edges of the MST in Step E) to the 
edges of an RST on the grid graph; [O(k)] ] 
[Theorem 1] : 
The time complexity of MAIN ALGORITHM is 
O(klogk), where k is the number of generating points.g 
(Proof) First, in Step A), Delaunay net is ob- 
tained by using the divide-and-conquer technique[2] 
with O(klogk). 
Next,in Step B), the check for each triangular 
Steiner point is done in each triangle of the Delau- 
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nay net. Since the number of edges in Delaunay net 
is at most 3k and from Lemma 3, the number of th 
triangular Steiner points in Delaunay net is at nok 
k-2. Hence,the time complexity of Step B) is O(k) 

In Step C), the modified Delauney net of Py 5 te 
obtained by the same way as Step A),where however all 
the triangular Steiner points are also considered to be 
virtual generating points. Hence, the time complexity 
of Step C) is O(k'logk’), that is to say, O(klogk). (k!. 
the number of elements of PU S and k’ = 2k ~ 2) 

In Step D), on the modified Delaunay net an MS? 
of PU S is obtained by [ALGORITHM 1] or [ALGo. 
RITHM 2] with O(klogk) by Lemma 4 or Lemma 5 re. 
spectively. 

In Step E),deletion of the triangular Steiner point 
with the degree two or less from MST(in Step D) is 
performed by deleting the point and the edges incident 
with it. Since the number of points on MST is O(k). 
In this case, the time complexity is O(k). 

Step F) depend on the number of points on MST. 
In this case, the time complexity is also O(k). 

Therefore [MAIN ALGORITHM] can be calculated 
in O(klogk). i 


[ALGORITHM 1] 
Starting Condition ) 
M := 1; /* Label Number */ 
Tt = ¢; 
/* T ; The subtree of an MST on U = PU S*/ 
7 A te 
/* TA; point set of subtree of an MST on the modified 
Delaunay net */ 
Start) Choose a; from P. 
TA := {a1}; 
U := U — {ai}; 
Label ai as M. 
Step 1) Choose a minimum weight edge 'E' among 
edges incident with the point M and any point u €U 
under [ CONDITION 1]. 
If there does’t exist ‘EB’ satisfying [CONDITION 1}; 
then go to Step 6. 
Step 2) Choose a minimum weight edge between any 
point belonged 7’A and the point 'Q’ under [CONDI- 
TION 2]. /* Q ; another end point of the edge 'E’ | 
Step 3) Connect the edge[Q, M] between Q and M, 
where M is the end point of the minimum edge se 
lected in Step 2) and is not Q. 
TA :=TA+{Q}; 
U := U - {9}; ; 
add the point Q and the edge [Q, M] to T; 
M:=M+4+1; 
Label Q as M. 
Step 4) If U = ¢, then stop. 
Step 5) Go to Step 1. 
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[CONDITION 1] 

In case of two or more edges with the minimum weight, 
we have to choose the edge with another end point 
which is the Steiner point and is also an element of 
U. Moreover, if there exist two or more edges satis- 
fying the above condition, we have to choose one of 
them smallest first with respect to abscissas of them. 
If there exist two or more edges satisfying the right 
above condition, we have to choose one of them small- 
est first with respect to ordinates of them. 
[CONDITION 2] 

In case of two or more edges with the minimum weight, 
we have to choose the edge with another end point 
which is the Steiner point and is also an element of 
TA. Moreover, if there exist two or more edges satis- 
fying the above condition, we have to choose another 
end point which is the Steiner point with the maximum 
degree. If there exist two or more edges satisfying the 


‘right above condition, we have to choose one of them 


as the same way as [CONDITION 1]. 

[SUB ALGORITHM-1] 

Step 6) J:=1; 

Step 7) M:=M-1; 

Choose a minimum weight edge ‘EF’ among edges in- 
cident with the point M and any point u € U under 
[CONDITION 1]. 

If there does’t exist 'E' satisfying [CONDITION 1], 
then go to Step 9. 

Step 8) M:=M4J; 

Go to Step 2. 

Step 9) J:=J+41; 

Go. to Step 7. i 

[Lemma 4] 

The complexity of »ALGORITHM 1” is O(klogk). 
q 
(Proof) For Step 1 and Step 2, the number of searching 
edges are both O(k) and the complexity of these steps 
is O(k) as well. 

In the case of Step 3, in order of having a loop that 
accounts for the total number of generating points and 
triangular Steiner points, (all of which are elements of 
set U(U = PUS)), the complexity of this step has to 
depend on the elements of U. 

In the case where the number of points is O(k), the 
time complexity of this step is O(k). Similarly, the 


‘time complexity of Steps 4,5, and 6 is O(k). 


The time complexity of Steps 6,7,8, and 9, according 
to the divide-and-conquer technique, is O(klogk). 

Accordingly, the entire algorithm( ALGORITHM 1 
) is O(klogk). I 
[ALGORITHM 2] 
Step 1) Sort all the points of U = PUS smallest 
first with respect to abscissas of points and if there are 
Many points with the same abscissa, sort them small- 


est first with respect to ordinates of them. [O(klogk)] 
Step 2) Label each point according to the point se- 
quence obtained in Step 1. [O(k)] 

Step 3) Sort all the edges of the modified Delaunay 
net smallest first with respect to weights of edges and if 
there are many edges with the same weight, sort them 
smallest first with respect to their end point labelings. 


[O(klogk)] 


‘Step 4) Label each edge according to the edge se- 


quence obtained in Step 3. [O(k)] 

Step 5) Make a cluster which is composed of two 
points selected successively from U in ascending label 
under [CONDITION] and label it one by one. 

If [CONDITION] is not satisfied, make a Culter con- 
sists of just one point with smaller label and proceed 
the clustering by initiating the other point with lager 
label. If it remains just one point to pair, this one 
point is the last cluster. [O(k)] 
Step 6) Connect the edge between two points in each 
cluster and delete all data about all the connected 
edges. [O(k)] 

Step 7) Connect the edge with the smallest label be- 
tween two clusters specified by this edge and make new 
cluster having the smaller label number of the above 
two old clusters and dump all the edge data except _ 
subtree edges within the new cluster. If it remains 
just one cluster with the label number ”1”, then stop. 
[0(k)] | ‘ 

[CONDITION] There exists an edge between two 
points on the modified Delaunay net. a 
[Lemma 5] 

The time complexity of ALGORITHM 2 is 
O(klogk), where k is the number of generating points. 


(Proof) We can use a sorting algorithm called heap- 
sort for Step 1 and Step3, then each time complexity 
is O(klogk). In Steps 2, 4, and 5, the labeling for sub- 
set of points or edges of net is proceeded, then each 
time complexity is O(k). Moreover, Step 6 and Step 7 
are treating points of net, then each time complexity 
is also O(k). Therefore we have conclusion. i 


5 EXAMPLES 


EXAMPLE 1(THE FIRST METHOD USING ALGORITHM 1) 


e Let us consider that six generating points are specified as shown in Fig.1- 
1. By Step A and Step B of MAIN ALGORITHM, we have Fig.1-1 with 
Delaunay Triangulation and four triangular Steiner points. 


By Step C, we have the modified Delaunay net with points | U |= 10 as 
shown in Fig.1-2. 


If we select p(1) as the starting point of Fig.1-2 and proceed Step 1 to 
Step 3 of ALGORITHM 1 for Fig.1-2, we get Fig.1-3, where the set T 
is subtree and p(6) is labeled as @for convenience. 


In Fig.1-3, p(5) is the candidate for © by Step 1. But, by Steps 2 and 
3, ©is selected as M and we get Fig.1-4. Moreover, we have Fig.1-5 by 
repeating Steps 1 to 3. In Fig.1-5, the point ®has no candidate, then 
we have Fig.1-6 by using SUB ALGORITHM-1. 
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e Repeating again Steps 1 to 3 for Fig.1-6, we have finally an MST on the 


modified Delaunay net as shown in Fig.1-7. 


e By applying Steps E and F of MAIN ALGORITHM to Fig.1-7, we get 
an RST with cost 11 as shown in Fig.1-8. This is the answer of the first 


proposed algorithm. 

6 CONCLUSION 

In this paper, first, we have proposed two new 
approximation algorithms with the time complexity 
O(klogk),k is the number of the given generating 
points, for RST problem, which both algorithms make 
use of the modified Delaunay net defined by the given 
generating points and the derived triangular Steiner 
points. In order to make successively a minimum span- 
ning subtree on the modified Delaunay net, we used 
doubly the neighborhood property [5] of each triangu- 
lar Steiner point [6] as well as each generating point for 
the first method, and we used Kruskal’s Algorithm-like 
procedure [8] for the second method. 

It is also natural to say the first method is the ex- 
tended version of Ref.[6], but, in Ref.[6], the full in- 
formation contained in the modified Delaunay net has 
not been used to make an MS subtree, while in the 
first method as well as the second one of this paper, 
it is used thoroughly to make an MS whole-tree, not 
subtree. We think that the proposed two algorithms 
are best among the known approximation methods 
(4],[5],[6] with respect to both the time complexity and 
the solution weight, because the local information to 


ified Delaunay net. Details about these facts will be 
shown somewhere. 
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